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We analyse a model of a one-dimensional electron gas interacting with an antiferromagnetic Heisen- 
berg spin-1/2 chain via the spin exchange interactions. Using a solution at a special limit, we 
characterize the gapless modes of the spin gap fixed point at weak coupling Jk <C Jh,Ef- we 
show that the only gapless pairing mode with divergent susceptibility is a composite odd-parity 
odd-frequency singlet pairing order parameter, while the ordinary BCS even-parity singlet pairing 
mode is incoherent. For 2-leg ladder systems, we note that it is possible to have a range of doping 
where the chemical potential cuts only the anti-bonding band while the bonding band remains half- 
filled. We propose that, in such a state, the 2-leg ladder is effectively realizing the one-dimensional 
Kondo-Heisenberg model. 



The one-dimensional Kondo-Heisenberg model (K-H 
model) describes an incommensurate one-dimensional 
electron gas (1DEG) interacting with a Heisenberg chain 
of spins- i via spin exchange interaction. For the K- 
H model, we show that the only gapless pairing mode 
with divergent susceptibility is a composite odd-parity 
odd-frequency singlet pairing order parameter, while the 
ordinary BCS even-parity singlet pairing mode is inco- 
herent. In addition, we find that the generalized Lut- 
tinger's theorem of Yamanaka et al. |Q, is satisfied only 
by the introduction of a new composite charge density 
wave (CDW). The composite CDW has power-law cor- 
relations with a " large- Fermi-sea" characteristics, while 
conventional CDW correlations decay exponentially. 

We discuss the significance of our results in several con- 
texts: First, our analysis sheds new light on the relations 
between previous treatments |2||j| of the one-dimensional 
Kondo-Heisenberg model. Second, we discuss the pos- 
sibility of effective realization of K-H model physics in 
2-leg ladder systems (This possibility was missed in all 
previous studies of doped 2-leg ladders g). Third, we 
criticize previous suggestions regarding the relevance of 
K-H model to "stripes theories" of high-Tc superconduc- 
tors. 

The core of our analysis is based on the derivation of a 
special solvable limit of the K-H model, and the meaning 
of such a solution within the renormalization group (RG) 
framework. The previous perturbative RG analysis of the 
K-H model Q has shown that the spin exchange interac- 
tion flows to some strong coupling fixed point suggesting 
the formation of a spin gap phase with enhanced pairing 
correlations. For particular value of parameters we ob- 
tain a well-controlled analytical solution which enables us 
to enumerate and characterize quantum numbers of all 
gapless modes. Gapless modes are properties of the fixed 
point. This means that the same gapless modes char- 
acterize all models which flow to the same fixed point. 
In particular, if there is only one fixed point to which 
all weak coupling K-H models flow, then our analysis is 
valid for all of them. 

The K-H model (0) consists of two inequivalent inter- 



acting chains; one is a one-dimensional electron gas (de- 
scribed by the Hamiltonian H 1DEG ||), and the other an 
antiferromagnetic Heisenberg chain of localized spins 1/2, 
{fj}. The chains interact via a spin exchange interaction 
with an antiferromagnetic coupling constant Jk > 0. 



H = H lu ^ +H Hels +H K (1) 

jjHeis =Jh Y j fj ■ fj+l , H K = J K Tj ■ S (Xj ) (2) 

3 i 

where s(xj) = ^a(xj)^-^)ff(xj) is the electron gas 
spin density operator at position Xj of the local spin 
fj of the Heisenberg chain. We focus on the low en- 
ergy and long-distance behavior of the electron's corre- 
lation functions by taking the continuum limit of the 
electron gas and linearizing the 1DEG dispersion rela- 
tion about the fermi points, ifcp, with corresponding 
right and left going electron fields, R a and L a ; tp a (x) — 
R a {x)e +ikFX + L a (x)e- tkFX . Where a =f, |. 

The 1DEG spin currents are decomposed into forward 
and back-scattering parts; 



s i x ) = ^l(Xj)^-^0(Xj) 



= [3 sR (x) +3 sL (x)] + n s (x) 



where 3 S R = hRta 



J S L = ^L+(j aa i L n > are the 



2 R>(j @ aa' R> 

ferromagnetic (q — 0) spin currents of right- and left 
moving electrons respectively, and 



n s (x) = 
where n# = R£ 



>n R (x)+e +i2kFX *n L (x) 



(4) 
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gered magnetization {q = 2Uf) components of the 1DEG. 
We work in the weak inter- chain coupling limit 



—R a ' are the stag- 



Jk <C Jn,Ep. 



(5) 



In this case, one is allowed to make further approxima- 
tion by taking the continuum limit also for the Heisenberg 
spin chain Q (such approximation is not valid in the op- 
posite limit Jk Jh, which is discussed elsewhere HH). 
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The local spin chain field is then also decomposed into 
the smooth (ferromagnetic) and staggered (antiferromag- 
netic) components; 



Tj = [3r (xj) + 3 T L (xj)] + (-I) n r ( Xj ) . 



(6) 



(Note: we will consistently use the subscripts "r, s" to 
distinguish the spin chain fields from the 1DEG fields). 

The effective Fermi wave numbers (in the sense of the 
generalized Luttinger's theorem for the 1DEG and 
the spin chain are 2kp and 2kp els = ir/b respectively 
(where b the distance between the local 

spins of the Heisenberg chain). It is assumed that the two 
systems are relatively incommensurate, and that 2kp is 
incommensurate with any underlying ionic lattice. In 
order to distinguish contributions coming from various 
interaction terms, we introduce distinct Kondo coupling 
coefficients for forward scattering (J/), back-scattering 
(Jb) and mixed interactions (J m ); 

H K = j f {r R + r L )-{y R + r L ) (?) 

+J m (-l) j n T - (J s fi + J£) 

+J b T(x j ) ■ [e' akF ^n sR (ay) + e +akFX ^n sL (x)] 

The back-scattering term J b and the mixed interaction 
J m are made irrelevant by the oscillatory factors e^ t2kFX: > 
and (— 1) J respectively. Therefore, at incommensurate 
filling in the weak coupling limit, the K-H Hamiltonian 
M reduces to 



H = H + J f J dx (3 R + J£) • (3% + 3D 



(8) 



where H Q = H 1DEG + H Hels . Due to the incommen- 
surate electron filling, after dropping terms which are 
irrelevant in the renormalization group (RG) sense, the 
spin and charge sectors decouple, H = j dx [Ti c + H S pi n }- 
The charge sector is described by a Gaussian model B 



Vc 
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K c n 2 c (x) + ±- (d x <f> c f 



(9) 



The subsequent analysis and manipulations deal only 
with the spin sector fields. The spin part of Hq can 
be written as the sum of two level k = 1 SU(2) 
Wess-Zumino-Novikov-Witten models. The correspond- 
ing Hamiltonian density is 



(10) 



where v T , v s are the spin wave velocities of the Heisenberg 
chain and 1DEG respectively (y T = ttJh/2). 

For Jf > Q, the non- interacting fixed point ([To]) is un- 
stable M , and the low-energy physics is governed by some 
"strong-coupling" fixed point. Nothing more can be de- 
duced from perturbative RG analysis, and the character 
of the strong coupling fixed point should be studied by 



means of non-perturbative methods. Sikkema- Affleck- 
White jl noticed that the relevant spin sector of the 
K-H Hamiltonian at incommensurate filling is equivalent 
to that of the 2-leg zigzag spin ladder. In turn, the zigzag 
ladder was shown to possess a spin gap by means of exact 
numerical simulations. 

We use the bosonized representation of the spin-i 



fermionic fields 



x) 



-i^Tr[e a (x)-<j} tJ (x)} _ Where 



R a (x) = F s=e -W^.W-<PAx)\. Where 0Jx) 

J^ oo dx'U a (x'), and [IL„(x'),</) a (x)] = -iS(x' - x), a =] 
The Klein factors, {F a ,F a i} = 5 a ^>, enforce proper 
anticommutation of fermions with different spin. As 
commonly done, we re-express the operators in terms 
of bosonic spin fields 4> s (x) — -^[<^f — ^M.]> anc ^ charge 
fields 4> c {x) = -4= [4>-\ + 4>i\, and correspondingly defined 
momenta H s and n c . Similarly, the spin chain fields are 
bosonized. In particular, the bosonized expression for 
the staggered magnetization is 



sin(v 27T# T ), — cos(v 2n8 T ), sin(v 2tt(/) t ) ) . (11 



In what follows we shall also need the bosonized expres- 
sion for the triplet superconducting order parameter: 



2 

~ e 



(12) 



\s/2^9 c (sin{s/2^6 s ), - cos(V2^6» s ), - sin(V27r0 s 



The model (g) is exactly solvable by the Bethe ansatz 
for general v s ,v r . However, for the purpose of cal- 
culating correlation functions, we take advantage of a 
special point in parameter space where the spin velocities 
are equal, i.e., Av s = v s — v T = 0. 

Using a transformation to composite spin fields; 9± = 
(0s i Qt) an d 4>± — -4^ (0s ± T ), the spin sector of 
the Kondo-Heisenberg Hamiltonian is simplified to the 
form; H = H§ z + AH ZZ + H_l 



H. 
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dx 



H: 



n: 
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2-kv, 



(d x <!>-\ 



{dx4>-Y 



AH zz = ^ij dx{n+n_ + {d x <j> + ) (d x <j>.)} 



H t = 



J 



dx cos 



■ ^V47T0 + ^ + cos {V^n(j)^j 



(13) 



(14) 



(15) 



Where v = ^(v s +v r ). Note that the J z part 
of the interaction has been completely absorbed into 
the kinetic energy part ( |l5| ) of the Hamiltonian in 
terms of the new fields, tfi±. Intuitively, a spin gap 
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can be established due to the j{ interaction term 
cos (\/47r#-) cos (\/47r</ , +) m (P"5[h where a self-consistent 
expectation value can be obtained for the composite fields 
(cos (V47T0-)) 7^ and (cos ( V^4>+) ) 7^ 0. In con- 
trast, (cos(V47T0_) cos(a/4^>-)) = (since e iV ** e - 

and e 4V ^ r ^- are respective disorder/order parameters) 
|j . Therefore, for determining spin gap physics of 
the fixed point, we rigorously need to keep only the 
cos (\/AkQ-^ cos (V*7r^+) interaction term. An addi- 
tional simplification (which we justify later on) is ob- 
tained if we neglect the velocity renormalization in (|l3| ) 
(i.e., equivalent to the anisotropic J z = limit). Thus, 
we obtain a decoupling of the spin sector into two com- 
muting sine-Gordon type Hamiltonians, 



A 



a — 1 n 



+ (-if ^ (d x ^) (d x e. t ) + cos 

4 2 (ira) 



(16) 



where, $i are new non-chiral fields combining the chiral 
components of <f> s and <j> r as follows: 



$1 = 



9i = 



V2 



V2 



; $2 = 



; e 2 = 



V2 



V2 



In the limit Av s = 0, the Hamiltonian ( |16[ ) is equivalent 
to the spin sector of the SU(2) Thirring model which is 
known to have an exponentially small gap |@), as antici- 
pated by the RG arguments . The spin gap fixed point 
is perturbatively stable with respect to all the interactions 
which were neglected for arriving at (jilq). In the ground 



(17) 



state \Z2n^>j = irn, where n is an integer. Thus 



cos(V27T$j) ) f- 0, ( sin(V27r$ i ) ) = 



and there is an additional discrete Z 2 x Z 2 symmetry cor- 
responding to the signs of (cos(v / 27r < I > :( )) which is sponta- 
neously broken in the ground state, and to the (<E>i,<I> 2 ) 
separation (the later is only an approximate symmetry 
which is broken by J z =/= coupling terms). 

A spin gaped one-dimensional system is expected 
to manifest enhanced pairing and charge density wave 
(CDW) correlations. Furthermore, the generalized Lut- 
tinger's theorem |Q mandates the existence of a gapless 
CDW mode at wave- vector 2k* F = 2Uf + ¥• As we shall 
see, these intuitive expectations are satisfied in a rather 
non-trivial manner. 

The rigidity of the composite bosonic fields, enforced 
in (|l7j), implies that the correlation function of any or- 
der parameter for which the spin part cannot be writ- 
ten purely in terms of cos(v / 27r ( I > j) is exponentially de- 
caying, i.e., is incoherent. In particular, the usual 
one-dimensional electron gas singlet charge-2e pairing 



L\R\ 



±(R\L ]+ R\L i)+ n.c. 

Instead, there are gapless modes of a composite nature: 
A composite odd-parity odd-frequency singlet 



and the 2k F CDW C dw 
are incoherent! 



Oc-sp {x) = —\Rl (ffa 2 ) af3 Ll 



(18) 



>s (V2tt$i) cos (V2tt$ 2 ) ) (-1) J 
and a composite CDW (a charge-0 spin-0 operator) 



Oc-CDW [x) = niDEG ■ T, 



(19) 

(2k F X + 7Tj) 



~ e +iV ^^ <Jcos (v^tt*i) cos (\/2^$ 2 ) \ e +l 

where it is the staggered component of the impurity spin 
chain, t^(— 1)'' n T , which is contributing the gapless 
modes with power-law correlations. 

The staggering factor (—1)"' in the corresponding corre- 
lation functions is effectively modulating the usual corre- 
lations by the reciprocal lattice vector ^ of the spin chain. 
As a result, the composite gapless modes are found at 
unusual finite momentum values: the composite singlet 
pairs with momentum ^ (and there is no k = singlet 
pairing with charge 2e), and the gapless composite CDW 
mode at momentum 2k* F = 2fcp+f (and not at 2kp of the 
bare 1DEG). The pure charge sector is not affected, as is 
evidenced by the fact that the gapless r\ — pairing mode 
(t]r = R\R^ and j]l — L^L^) remains at momentum 2fcj?. 
The gapless modes are inter-related by the commutation 



relation 



Oc- 



CDW 



Vr-Vl 



V2 



20 c -sp- An extended dis- 
cussion of these order parameters can be found in |6| . 

To summarize, we developed a solution of the one di- 
mensional Kondo-Heisenberg model at weak exchange 
coupling Jk <C Jh,Ef-i for a special value of parame- 
ters v s — v T . We were able to explicitly demonstrate the 
spin gap and characterized the gapless modes properties 
of the fixed point previously alluded to by perturbative 
RG arguments. 

The spin wave velocities v s and v T in general are not 
identical. What then are the limits of validity (and hence 
the significance) of our solution? First, the spin gap guar- 
anties that our fixed point solution is perturbatively sta- 
ble for Av s — v s — v T 7^ 0. Moreover, unless there is 
a phase transition driven by large Av s anisotropy (to a 
yet another unknown phase), our results (gapless modes 
identification) in the limit Av s — are universal for all 
Av s so long as the weak coupling condition Jk <C f T , v s 
is maintained. The same argument applies to all other 
approximations we undertook for arriving at (|l6|). 

The small Jk condition manifests itself through the co- 
efficient ^1 — °f t nc (<9x0-) 2 term in the Hamilto- 
nian (|l3|). It indicates that something may indeed break 
down when J I > 2ttv s ||. Indeed, we remark in this 
context that J/ = 2ttv s is the so called " Toulouse point" 
value on which we further comment below. Therefore, 



3 



the perturbative RG flows to "strong coupling" found in 
should be understood as being only to some interme- 
diate coupling fixed point <// < 2irv s with a finite basin 
of attraction (the true strong Jr- coupling limit of the 
K-H model is gapless JhJ). 

We now discuss the significance of our results in var- 
ious context: In previous paper ||, a spin gap phase of 
the Kondo-Heisenberg model (0) was found in another 
region of the parameters space (Jh <C Jk ~ Ep); the so 
called "Toulouse limit" solution. That solution has the 
same composite order parameters as the present one, but 
in addition it also possesses gapless modes of the conven- 
tional CDW and even-parity singlet pairing A order pa- 
rameters. Hence, we conclude that the "Toulouse point" 
H and the " zigzag ladder limit" || spin gap phases of the 
Kondo-Heisenberg model are distinct. An elaborate com- 
parison and implications for the general phase diagram of 
the Kondo-Heisenberg model is presented elsewhere . 

The two-leg ladder system consists of two equivalent 
Hubbard model chains which are coupled by single par- 
ticle hopping with amplitude fj_. As is well known, the 
effect of t± interaction can be treated exactly by intro- 
ducing " bonding" and " antibonding" bands described by 
fermion fields %j) A B = -j= (tp 1 ± i/> 2 ), where ip 1 2 (%) are 
fermion fields on the legs {1,2} of the ladder. The bond- 
ing and antibonding bands are inequivalent. In particu- 
lar, there is a chemical potential difference (ia — (J>b ~ tx- 
Consequently there could exists a range of doping for 
which, depending on model parameters, holes may enter 
only into the antibonding band (which becomes gapless) , 
while the bonding band remains half-filled and retains a 
Mott-Hubbard gap. In that range of doping, the Fermi 
energy lies in the gap of the bonding band and cuts only 
the antibonding band. As far as the low energy physics of 
such a state is concerned, the half-filled band is equivalent 
to a Heisenberg chain of localized spins 1/2, {fj}, with 
effective antifcrromagnetic coupling Jh- On the other 
hand, the gapless band represents a one-dimensional elec- 
tron gas (described by the Hamiltonian H 1DEG ||), with 
an incommensurate Fermi momentum kp. Hence, the 
only relevant interaction between the two bands is spin 
exchange interaction, Jk > 0. We conclude that in this 
particular doping range, the low energy physics of the 
2-leg ladder is effectively captured by the K-H model 
(|l|). The conditions under which the above scenario is 
realized require an elaboration beyond the scope of the 
this paper. The exact dependence of model parameters 
{Jh, Jk} on the original ladder parameters {t, U, V, t±} 
is unimportant since, as noted in the introduction, our 
analysis addresses general fixed point properties. Con- 
sequently, we suggest that in the general phase diagram 
of doped ladder 0], the spin-gap region (labeled C1S0 in 
Pj) should be divided in two: A certain low doping re- 
gion with only odd-«7 pairing, and a higher doping region 
with conventional pair states as discussed for example in 
Consequently, the same would be true for the puta- 
tive superconducting state in a system of coupled 2-leg 



ladders. 

A model of an incommensurate 1DEG coupled with a 
ladder environment [jllj was proposed in the context of 
stripe phases in HTc cuprates. In the case of a gapless 
spin ladder environment, if only spin exchange interac- 
tions are considered jljj one arrives at the effective model 
given by ([!]). Combining our analysis (where only odd- 
lu pairing is coherent) with the experimental observation 
that superconductivity in high-T c cuprates is due to d- 
wave BCS paired electron, we conclude that, within a 
stripe state scenario spin exchange interactions are 
ruled out as a possible source of the spin gap in high-T c 
cuprates. While unrelated to cuprates, our solution does 
indicates the possibility of making pure odd-time com- 
posite pairing superconductors by constructing 2D or 3D 
weakly coupled arrays of the ID chain model (Q) . 
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